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We propose a fully quantum-mechanical method of treating four-body nuclear breakup processes in scatter- 
ing of a projectile consisting of three constituents, by extending the continuum-discretized coupled-channels 
method. The three-body continuum states of the projectile are discretized by diagonalizing the internal Hamil- 
tonian of the projectile with the Gaussian basis functions. For ®He-l-^^C scattering at 18 and 229.8 MeV, the 
validity of the method is tested by convergence of the elastic and breakup cross sections with respect to increas- 
ing the number of the basis functions. Effects of the four-body breakup and the Borromean structure of ''He on 
the elastic and total reaction cross sections are discussed. 

PACS numbers: 21.45.+V, 21.60.Gx, 24.10.Eq, 25.60.-t 



The study on neutron-halo nuclei has become one of the 
central subjects in the unstable nuclear physics since the dis- 
covery of such nuclei 1 1]. In scattering of a two-neutron-halo 
nucleus such as ^He and '^'^Li, the projectile easily breaks up 
into its three constituents (n-i-rt-i-core), indicating that the scat- 
tering should be described as a four-body (n-i-n-i-core-i-target) 
reaction. Then an accurate theory for treating such a four- 
body breakup is highly desirable. 

So far the eikonal and adiabatic calculations were pro- 
posed and applied to ^He and '^'^Li scattering around 50 
MeV/nucleon Since these calculations are based 

on semi-classical approaches, they work well at higher inci- 
dent energies. In fact, the elastic cross section of ^He-i-'^^C 
scattering at 229.8 MeV has recently been measured |6] and 
successfully analyzed by the eikonal calculation with the six- 
nucleon wave function of ^He |7]. However, these approaches 
seem not to be applicable for low-energy scattering such as 
i2C('^He,'^He)i2C at 3 MeV/nucleon Hi measured very re- 
cently. 

In this rapid communication, we present a fully quantum- 
mechanical method of treating four-body nuclear breakup. 
The method is constructed by extending the continuum- 
discretized coupled-channels method (CDCC) [£] that treats 
three-body breakup processes in scattering of the two-body 
projectile. In CDCC, the total scattering wave function is ex- 
panded in terms of bound and continuum states of the projec- 
tile. The continuum states are classified by the linear {k) and 
angular momenta, and they are truncated by setting an upper 
limit to each quantum number. The fc-continuum is then di- 
vided into small bins and the continuum states in each bin are 
averaged into a single state. This procedure of discretization is 
called the average (Av) method. The S'-matrix elements calcu- 
lated with CDCC converge as the modelspace is extended |9||]. 
The converged CDCC solution is the unperturbed solution of 
the distorted Faddeev equations, and corrections to the solu- 
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tion are negligible within the region of space in which the re- 
action takes place lHoll . 

Also for four-body breakup processes in scattering of the 
three-body projectile, CDCC has to prepare three-body bound 
and discretized-continuum states of the projectile. Because of 
the difficulty of preparing all the three-body states with the 
Av method, CDCC so far analyzed ^He scattering within a 
limited model in which a two-neutron pair is treated as a single 
particle, di-neutron (^n) lITlll . However, the accuracy of the 
di-neutron model has not been confirmed yet, because of the 
absence of fully quantum-mechanical method of treating four- 
body breakup. 

In our previous work ifl^l on three-body breakup in scat- 
tering of the two-body projectile, we proposed a new method 
of discretization, called the pseudo-state (PS) method. In the 
method, continuum states of the projectile are replaced by 
discrete pseudo-states obtained by diagonalizing the internal 
Hamiltonian of the projectile in a space spanned by the L^- 
type Gaussian basis functions. The CDCC solution calculated 
by the PS method agrees with that by the Av method which 
can be regarded as the exact solution. Thus, a reasonable num- 
ber of the Gaussian basis functions can form an approximate 
complete set in a finite configuration space being important 
for three-body breakup processes. It is very likely that the ap- 
proximate completeness persists also in the case of four-body 
breakup processes. Actually, as shown latter, we can see clear 
convergence of calculated elastic and breakup cross sections 
with respect to increasing the number of the Gaussian basis 
functions. It should be noted that the Gaussian basis functions 
are widely used to solve bound-state problems of few-body 
systems LL2J, since the use of the basis functions reduces nu- 
merical works much. Thus, the four-body breakup processes 
can be analyzed properly by CDCC with the PS method. We 
refer to this new method as four-body CDCC and the usual 
CDCC for three-body breakup as three-body CDCC. 

The first application of four-body CDCC thus designed is 
made for '^He-n^^C scattering at 18 and 229.8 MeV, where the 
projectile has the BoiTomean structure and then easily breaks 
up into two nucleons and '*He. In these scattering processes. 
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the incident energies £'in are much higher than the Coulomb 
barrier energy (~ 3 MeV), so only nuclear breakup processes 
become significant. We thus concentrate our application on 
nuclear breakup. The calculated elastic cross sections well 
reproduce experimental data at both Em- Moreover, effects 
of the four-body breakup and the Borromean structure of ^He 
on the elastic and total reaction cross sections are discussed in 
the case of Em =18 MeV. 

We assume that ^He+^^C scattering is described as a four- 
body system, n+n+'^He+^^C. Then, the Schrodinger equation 
can be written as 



+ Vc{R) + He~E 



vi'(e,R) = o, (1) 



where R and ^ are, respectively, the coordinate of the center- 
of-mass of ^He relative to ^^C and the internal coordinates 
of ^He; Kji is the kinetic energy associated with R. Here, 
He is the internal Hamiltonian of ^He-projectile, and E is the 
sum of -Ein and the ground state energy of ^He. The Vij rep- 
resent two-body nuclear interactions working between ^He- 
projectile (P) and ^^C-target (T). Meanwhile, the Coulomb 
potential Vq is treated approximately as a function of R only, 
i.e., we neglect Coulomb breakup processes. 
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FIG. 1: Jacobian coordinates of three rearrangement channels (c 
1-3) adopted for the n+n+'^He model of ''He structure. 



The four-body wave function vj/'^J^^^ where J is the total 
angular momentum of the four-body system and M is its pro- 
jection on the z-axis, is expanded in terms of a finite number 
of the internal wave functions $^ of ^He-projectile: 
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where y^f^^ - [$„/(0 ® i'^YL{nn)]jM- The 7 stands for 
the set of (n, /, m), where / is the total spin of ^He and m is 
its projection on the z-axis, and n stands for the n th eigenstate 
with positive energy. The ground state of ^He, which is the 
only bound state of ^He, is denoted by 70 = (0, /o, toq)- The 
satisfies Hq^^ = e„/$7 and the expansion-coefficient 
X'ni L i'^ '^l- Q represents the relative motion between the 
projectile and the target; L is the orbital angular momentum 
regarding R. The relative momentum Pni is determined by 
the conservation of the total energy: E — P^j/2ii + e„/, 
with fi the reduced mass between the projectile and the tar- 
get. Multiplying Eq. Q by ^/ from left, one can ob- 
tain a set of coupled differential equations for Xni L' called 



CDCC equation; it should be noted that the CDCC equation 
for the four-body system is formally equal to that for the three- 
body system. Solving the CDCC equation under the appropri- 
ate asymptotic boundary condition |9, 14], we can obtain the 
elastic and discrete breakup 5-matrix elements. Details of the 
formalism of CDCC are shown in Ref . 0| . 

In the Gaussian expansion method (GEM) fl3\, each is 
given by 



(3) 



where c denotes the set of Jacobian coordinates shown in 
Fig.n We here take the angular-momentum coupling scheme 
as I = A + S, where A and S are the total orbital-angular- 
momentum and the total intrinsic-spin of ^He, respectively. 

(c) 

Then has the following form: 



A.S ^ ^ 

(4) 

where rii/2 is the spin wave function of each valence neutron 
(ill or n2) and ''He has been treated as an inert core with the 
(Os)** internal configuration, (p''°'\ The definition of (y^, r^) 

is given in Fig.H The amplitude-function ^^aa/a' '^^^^ -^^A 
the projection of A on the z-axis, is expanded in terms of the 
Gaussian basis functions: 



^ilx.(yc,re) = E E E <UA2/^^e-(-/^-)^ 
xe i=i j=i 



where A {t) is the angular momentum regarding yc (rc). The 
Gaussian range parameters are taken to be a geometric pro- 
gression: 



(6) 
(7) 



The $^ is antisymmetrized for the exchange between n\ and 



n2\ we then have A, 



(2) 



must be 1 for c = 3. Meanwhile, the exchange between 
each valence neutron and each nucleon in "^He is treated ap- 
proximately by the orthogonality condition model |15]. The 
eigenenergies e„/ of ^He and the corresponding expansion- 
coefficients A^ ^^jgj^ are determined by diagonahzing iJg 

El. 

In the four-body CDCC calculation shown below, we take 
= 0+ and 2+ states for "^He; tt is the parity of ''He. We 



TABLE I: The maximum internal angular momenta and the Gaussian 
range parameters for each Jacobian coordinate. 
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3 0+ 
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0.1 


10.0 
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10.0 


1,2 0+ 
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0.5 
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0.5 
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3 2+ 
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0.5 


8.0 


1,2 2+ 
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0.5 
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0.5 


8.0 
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TABLE II: The number of the Gaussian basis functions, imaif ' and jmax^'' for set I, II and III. The corresponding number of the eigenstates of 
He, TVjJiaxi ™d the number of channels included in the CDCC equation, N^^^, are also shown (see the text for the details). 
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set I 


8 
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6 6 


204 


6 6 


288 


17 
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28 


39 


set II 


10 


8 


8 8 


352 


8 8 


512 


25 


32 


44 


64 


set III 


12 


10 


10 10 


540 


10 10 


800 


32 


42 


60 


85 



show in Table m the maximum values of the internal angular 
momenta, Amax, (max and Amax, and the Gaussian range pa- 
rameters, yi, ymax, ri and fmax, used in the calculation of 
It should be noted that most of them depend on and 
c, while in Eqs. (|5}-0 the dependence has been omitted for 
simplicity. In order to demonstrate the convergence of the 
four-body CDCC solution with respect to increasing the num- 
ber of the Gaussian basis functions, we prepare three sets of 
the basis functions, i.e., sets I, II and III. Each set is speci- 
fied by iinax and jmax''; again, the and (c)-dependence 
of them has been omitted in Eqs. One can calculate 

the total number of the eigenenergies of Hq, AAj^ax' t>y using 
Eqs. (|3}-(0 and the input parameters shown in Table |I] The 
values of ilaax , juiaK and N^^^y. for each set are shown in 
Tableim In the actual CDCC calculation for '^HeH-^^C scatter- 
ing at 18 MeV (229.8 MeV), high-lying states with e„7 > 12 
MeV (e„7 > 25 MeV) are found to give no effect on the elas- 
tic and breakup S'-matrix elements. Thus, the effective num- 
ber of the eigenstates of ^He, N^^^, is reduced much for each 
of sets I-III, as shown in Tablellll 

As for the coupling potentials in the CDCC equation, we 
adopt the double folding model 1 18] as follows: 

U^{^{R) = {NR + iNi)V^{^iR), (8) 

xvNN{E,p,s)drTdrpdnR, (9) 

where rp (rx) is the coordinate of a nucleon in the projectile 
(target) relative to the center-of-mass of the particle, and s = 
R + rx — rp . The quantum number ( represents n, I and L 
together, and the elastic channel, which has the incident wave, 
is denoted by Co = (0, Iq, Lq). The ground state density of 

12c, pZ^vt) ^ {^Z.\E]ii SivT - r,)lC2)^ where 
is the wave function of ^^C in the ground state, is calculated 
by the microscopic 3 a cluster model |19]. We in this study 
define the transition densities of ^He, p^f^,^"^*^, as 

As for the nucleon-nucleon effective interaction wnn, we use 
the realistic energy- and density-dependent M3Y (DDM3Y) 
interaction i20ll . Since the DDM3Y interaction is real, 
V^,^{R) has no imaginary part. Thus, we have multiplied 
V^,j~{R) by a complex factor Nn + iNj. In the present anal- 
ysis, we fix Njf = 1 and optimize Nj to fit experimental data 



for elastic scattering. It should be noted that in three-body 
CDCC calculation made before for ^Li scattering on various 
target nuclei L2ij.22], the prescription above was successful 
in reproducing experimental data. 

The convergence of the four-body CDCC solution is tested 
for ^He-H^^C scattering at 18 MeV. Figure|2]shows the energy- 
integrated breakup cross section, i.e., the sum of the cross sec- 
tions to all breakup channels, calculated with sets I-III. The 
results of set II and set III are in good agreement with each 
other, but the result of set I is somewhat different from them. 
Meanwhile, as for the elastic cross section shown in Fig. |3] 
the three sets give the same cross section shown by the solid 
line. Thus, the four-body CDCC solution converges with set 
II. Furthermore, we have confirmed that similar convergence 
is also seen with respect to extending ^max and f,nax- The 
optimum value of Nj determined from the measured elastic 
cross section is 0.5 at Ein — 18 MeV, which is the same as that 
for ^Li scattering at various Ein ll2ll l22il . It should be noted 
that all calculations shown in Figs.|2]and|3]use the same value 
of Ni. Also for ''HeH-i^C scattering at 229.8 MeV, we can 
see similar convergence of the elastic and energy-integrated 
breakup cross sections with respect to extending the mod- 
elspace. Comparison between the calculated and measured 
elastic cross sections is shown in Fig. 0] In this case the op- 
timum value of N] is 0.3. In Figs. |3l and |4] the dotted lines 
represent the elastic cross sections due to the single-channel 
calculation. Then, the difference between the solid and dotted 
lines shows the effect of the four-body breakup on the elastic 
cross section. For both Ei^, the effect is sizable, properties of 
which are discussed later. 
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FIG. 2: Angular distribution of the energy-integrated breakup cross 
section for ^He+^^C scattering at 18 MeV. The dotted, dashed and 
solid lines are the results of the four-body CDCC calculation with 
set I, II and III, respectively, of the Gaussian basis functions. 
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FIG. 3: Angular distribution of the elastic differential cross section 
for "^He+^^C scattering at 18 MeV. The solid and dotted lines show 
the results with and without breakup effects, respectively. The exper- 
imental data are taken from Ref. 
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FIG. 4: The same as in Fig.|3lbut for ^He+^^C scattering at 229.8 
MeV. The experimental data are taken from Ref. Q|. 

Recently, it was reported in Ref. fllll that the total reac- 
tion cross section for ^He-H^^^Bi is much larger than that for 
^LiH-^*'^Bi at similar energies relative to the Coulomb barrier 
energies because of the large El excitation strength of ^He 
to the continuum. Meanwhile, for ^He-n^^C scattering at 18 
MeV, the El excitation of ^He is negligible because Ei^ is 
much higher than the Coulomb barrier energy (about 3 MeV). 
As shown in Fig.|5] however, we find that 10% enhancement 
of the total reaction cross section is still left. The open cir- 
cles represent the total reaction cross sections for ^He-n^^C 
at 18 and 229.8 MeV calculated by four-body CDCC, while 
the filled circles show those for ^Li-n^^C in the energy range 
20-318 MeV calculated by three-body CDCC |21,'2^, where 
the microscopic d+^He model is assumed for ^Li structure. 
As mentioned above, the resulting optimum Nj value for 
^Li-H^^C scattering is about 0.5, i.e., almost independent of 

Ein- 

In order to investigate the origin of the 10% enhancement, 
we perform the three-body CDCC calculation by assuming 
the di-neutron model for ^He structure; in the model, the 
di-neutron density is assumed to be the same as that of the 
deuteron, and then the resulting ^He density is close to the 



FIG. 5: The incident-energy dependence of the total reaction cross 
section for scattering of ^He and ''Li on ^^C. The open circles show 
the results for ^Ue+^^C scattering at 18 and 229.8 MeV calculated by 
four-body CDCC, while the filled circles represent those for "Li-l-^^C 
scattering at several energies calculated by three-body CDCC based 
on the d-l-*He model for ®Li structure. The open triangle is the result 
for ^He-l-^^C at 18 MeV calculated by thi-ee-body CDCC with the di- 
neutron model for "He structure. The filled triangle is based on the 
same calculation as the open triangle, except that the Coulomb po- 
tential between ®He and ^^C is replaced artificially by that between 
''Li and ^^c. 



^Li one. The result of this calculation is shown by the open 
triangle in Fig. |5l The difference between the open trian- 
gle and the open circle at 18 MeV is due to the Borromean 
structure of ^He, which is referred to as the Borromean ef- 
fect. The effect dominates about half the 10% enhancement. 
The rest of the enhancement is mainly due to the difference of 
the Coulomb barrier energies between ^He-n^^C and ^Li-n^^C. 
Actually, when the Coulomb potential for ^He-n^^C is re- 
placed artificially by that for '^LiH-^^C, the CDCC calculation 
based on the di-neutron model (the filled triangle) gives the 
total reaction cross section close to the filled circle at 20 MeV. 

As for ^He-H^^C scattering at 229.8 MeV, we have con- 
firmed through the same analysis that the Borromean effect 
becomes negligible as well as the effect of the difference of 
the Coulomb bamer between '^HeH-^^C and •^Li-n^^C. This 
suggests no enhancement theoretically. Nevertheless, Fig. |5l 
shows that the total reaction cross section for ^He-n^^C is even 
smaller than that for *'LiH-^^C at the similar energy. This cu- 
rious behavior is due to the fact that Nj — 0.3 for ^He-n^^C 
while N] = 0.5 for ^Li-n^^C at this high energy. In fact, the 
total reaction cross section is enhanced by changing Nj from 
0.3 to 0.5 in four-body CDCC calculation for ^He+^^C, and 
the resulting cross section almost reproduces the correspond- 
ing one for ^Li+^^C. The origin of the small Nj value for the 
^He scattering is not clear at this moment, so more systematic 
experimental data are highly desirable for ^He scattering. 

Finally, we calculate the dynamical polarization (DP) po- 
tential induced by the four-body breakup processes, in order 
to understand effects of the processes on the elastic scattering. 
The DP potential C/^p is given by 

[/f^p(i?) = C//q(i?)-C//„^„(i?), (11) 
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FIG. 6: The dynamical polarization potential for ^He+^'^C scattering 
at 18 MeV with the grazing angular momentum Jgr = 10. The solid 
and dotted lines, respectively, represent the real and imaginary parts 
of the DP potential calculated by four-body CDCC. The dashed and 
dot-dashed lines correspond to those of three-body CDCC with the 
di-neutron model for ®He structure. The dot-dot-dashed line repre- 
sents the imaginary part of the double-folded potential U^^^^ . 

Ui^{R) + Vc{R)-E,, 

XXo^/o,Lo(^0/o,i?) = 0,(12) 



where Xoig Lq scattering wave-function in the elastic- 

channel calculated with CDCC, and [7^^ defined by Eq. il2\ 
is so-called the equivalent local potential. The detailed def- 
inition of the DP potential is shown in Ref. ll2lll . Figure |6l 
shows the DP potential for the ^He-n^^C scattering at 18 MeV 
with the total grazing angular momentum Jg^ = 10. The "in- 
sensitive region" of R shown in the figure is defined with the 
condition that |xo/o Lo^^oIoj ^)I is less than 5% of its max- 
imum value in the asymptotic region. The DP potential is 
almost independent of J around Jgi in the peripheral region. 

In Fig. |6lthe real part V^p^"* (V^p"^^) and the imaginary part 

^DP ^ (W^DP ^) of C^DP calculated by four-body (three-body) 
CDCC are, respectively, shown by the solid (dashed) and the 



dotted (dot-dashed) lines. Both of V^p and V^p^' are repul- 
sive and have almost the same strength which is about 30% of 



-.7(3) 



the real part of /7^^^^^,. The W^p ^ is about 20% of the imag- 
inary part of t^^Q^j, (dot-dot-dashed line), while Wjjjp'^'' oscil- 

lates with R, so the net effect of Wj~,p is negligibly small. 
Thus, one sees that inclusion of the four-body breakup pro- 
cesses, i.e., beyond the three-body breakup, makes the real 
part of t7|/^„ slightly shallow and the imaginary one deep. In 
particular, the latter effect is important and can be assumed to 
come from the Borromean structure of ^He. This is consistent 
with the fact that the total reaction cross section is enhanced 
by the Borromean structure. 

In conclusion, a fully quantum-mechanical method of treat- 
ing four-body nuclear breakup is presented by extending 
the continuum-discretized coupled-channels method. The 
method called four-body CDCC is appUed to ^He-n^^C scat- 
tering at 18 and 229.8 MeV in which ^He easily breaks up 
into two neutrons and ^He. In four-body CDCC, three-body 
continuum states of ^He are discretized by diagonalizing the 
internal Hamiltonian of ^He with the Gaussian basis func- 
tions. The validity of four-body CDCC is confirmed by clear 
convergence of the calculated elastic and energy-integrated 
breakup cross sections with respect to increasing the number 
of the Gaussian basis functions. We can say from the conver- 
gence that the Gaussian basis functions form an approximate 
complete set in a finite configuration space being important 
for four-body nuclear breakup processes. Furthermore, we 
find a 10% enhancement of the total reaction cross section of 
•"He-H^^C at 18 MeV relative to that of '^Li-n^^C at the similar 
energy. Half of the 10% enhancement is due to the Borromean 
structure of ®He. For the elastic scattering, the four-body 
breakup processes make, in particular, the imaginary part of 
the double-folded potential deep, which is originated in the 
Borromean structure of ^He. In the present analysis four-body 
Coulomb breakup is neglected. However, it would be possible 
to treat the Coulomb breakup within the present framework, if 
the complex-range Gaussian basis functions are taken I23ll . 
Further work along this line is highly expected. 
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